Sound/elastic-wave systems with fermionic quantum-like behavior for phase-based encoding of information

ABSTRACT

Embodiments for sound and elastic-wave systems with fermionic quantum-like behavior and in particular to an approach for the phase-based encoding of information are disclosed.

FIELD

The present disclosure generally relates to sound and elastic-wave systems with fermionic quantum-like behaviour and in particular to an approach for the phase-based encoding of information.

BACKGROUND

Wave functions possess a spinor and a spatio-temporal part. The spinor part leads to a non-conventional torsional topology of the wave function. In the long-wavelength limit, field theoretical methods can be used to demonstrate that rotational phonons can exhibit fermion-like behavior.

The understanding of elastic waves has been largely driven by the simple paradigm of the plane wave and its periodic counterpart (the Bloch wave) in periodic media. Significant progress has been achieved in unraveling the behavior of elastic waves in phononic crystals and acoustic metamaterials. However, most of this progress has been achieved through exploration of the geometrical complexity of these media. Inspired by the discovery of topological insulators for which the electronic wave function is supported in momentum space by manifolds with non-conventional topologies, recent studies have shown the possibility of achieving electromagnetic waves and acoustic waves with non-conventional topologies.

It is with these observations in mind, among others, that various aspects of the present disclosure were conceived and developed.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 depicts a schematic illustration of a discrete micromechanics model that supports rotational waves.

FIG. 2 depicts a schematic illustration of the torsional topology of the wave function for rotational modes in k space.

FIG. 3 depicts a mass-spring system isomorphic to a block-spring system supporting rotational waves only.

FIG. 4 depicts a schematic illustration of a device composite system comprising a segment of a mass-spring system with fermion-like behavior sandwiched between two semi-infinite harmonic chains.

FIG. 5 depicts the coefficient of transmission of a plane wave launched in the right semi-infinite chain through the composite system of FIG. 4, with an inset depicting a magnified region around the bottom of the positive branch of the band structure.

DETAILED DESCRIPTION

The present application investigates the properties of wave functions in phononic structures supporting rotational waves which can exhibit wave functions that take the form of spinors. The non-conventional torsional topology of the wave function of elastic waves with spinor characteristics leads to a constraint on the wave function reminiscent of the Pauli-exclusion principle. This type of constraint introduces the notion of fermion-like behavior of elastic waves. The presently discussed topological constraints on elastic waves can improve phononic systems, permitting unique, robust designs and new device functionalities by providing immunity to performance degradation caused by imperfections or information coding and processing in the phase of the waves.

Previous investigations have shed light on various properties a 2D phononic crystal (PC) constituted of stiff polymer inclusions in a soft elastomer matrix. The 2D PC composed of a square array of polystyrene (PS) inclusions in a polydimethylsiloxane (PDMS) elastomer matrix has been shown to be capable of supporting rotational waves. Of particular interest are modes where the PS inclusions and the region of the matrix separated by the inclusions rotate out of phase but also in phase. Following Peng et al., who demonstrated that a 1D lumped mass model could be used to describe rotational modes in a 2D PC, a 1D mass spring phononic structure that could also support rotational waves was introduced. This 1D model can be used to reproduce the dispersion relations of the 2D system in a certain range, and it can subsequently be shown within Dirac's formalism that the wave function for rotational waves possesses a spinor-part and a spatio-temporal part.

The present application discusses the geometric topology of the spinor-part of the wave function, and shows that the spinor-part imparts a non-conventional topology to the wave function. In particular, the spinor part of the elastic wave function is supported in momentum space by a manifold with torsional topology. Furthermore, in the long-wavelength limit of the 1D model, the concept of phononic structures that support rotational waves and the concept of a Dirac representation of elastic waves within a field theoretical framework are integrated with each other. This leads to the observation of fermion-like behavior resulting from the spinor part of the wave function. In particular, it is found that within the context of a “second quantization” of the rotational waves, the properties of their associated spinor lead to anticommutation rules for the creation and annihilation operators.

After showing the possibility of fermion-like constraints on the elastic waves, a second question related to the more practical aspect of this behavior is also posed: what are the implications of fermion-like phonon behavior? A partial answer to that question is provided by introducing the concept of phase control of phonons through the specificity of their fermion character. To that effect, the concept of a phononic structure-based ψ(phase)-bit is introduced, and it is demonstrated that one can operate on the ψ-bit through the structure's physical parameters to transform the spinor state of the wave function. More specifically, one can operate on and measure the spinor state of phonons replicating known operations that are more commonly executed on fermions states like spins.

The 1D discrete mass-spring model that supports rotational modes forms the foundation of a field theoretical representation of rotational modes in phononic structures. The wave equation associated with that model is shown to be isomorphic to the Klein Gordon equation, and that equation is therefore factored using Dirac formalism to obtain the wave functions and investigate their non-conventional topology in momentum (k) space. Quantum field theoretic methodologies are used to analyze the properties of the wave functions in the long-wavelength limit and find constraints on the wave function reminiscent of fermion-like behavior. An application of these findings to the encoding and processing of information in the wave function is disclosed.

Nonconventional Topology of Rotational Waves in Phononic Crystals

One-Dimensional Discrete Micromechanics Model

The various embodiments makes use of a 1D mass-spring model that is known to reproduce the dispersion relations of a 2D phonon crystal that supports rotational waves. This model is based on a discrete linear one-dimensional micromechanics model that includes longitudinal, shear and rotational degrees of freedom. This 1D discrete lattice model comprises an infinite chain of square block elements connected with multiple harmonic springs. Each element in the model is considered to have two translational degrees of freedom (displacement in the x and y directions) and one rotational degree of freedom (rotation about an axis perpendicular to the xy-plane).

FIG. 1 depicts the repeatable unit cells for a monoblock lattice model 100 with periodicity (h). More specifically, FIG. 1 is a schematic illustration of a discrete micromechanics model that supports rotational waves. Three different harmonic springs (spring constants k₀, k₁, and k₂) connect different parts of the block elements 101. The element in FIG. 1 has mass (m) and moment of inertia (I). The block constituting the nth unit cell has x-displacement (u_(n)), y-displacement (v_(n)) and rotation component (φ_(n)). u_(n) and v_(n) represent displacements associated with longitudinal and transverse vibrations, respectively. The potential energy associated with the elastic connections 102 of elements (n) and (n+1) in the monoblock chain is written as follows:

$\begin{matrix} {{E_{n,{n + 1}} = {{\frac{1}{2}{K_{0}\left( {u_{n + 1} - u_{n}} \right)}^{2}} + {\frac{1}{2}{K_{1}\left\lbrack {\left( {v_{n + 1} - v_{n}} \right) + {\frac{h}{2}\left( {\varphi_{n + 1} + \varphi_{n}} \right)}} \right\rbrack}^{2}} + {\frac{1}{2}{K_{2}\left( {\varphi_{n + 1} - \varphi_{n}} \right)}^{2}}}}\mspace{79mu}{where}\mspace{79mu}{{K_{0} = \left( {\frac{k_{0}}{h^{2}} + \frac{2k_{1}}{I^{2}} + \frac{2k_{2}I^{2}}{I_{e}^{4}}} \right)},{K_{1} = \left( \frac{2{k_{2}\left( {2a} \right)}^{2}}{I_{e}^{4}} \right)},\mspace{79mu}{K_{2} = \left( \frac{2a^{2}k_{1}}{I^{2}} \right)},{I = {h - \left( {2a} \right)}},{I^{d} = {\sqrt{\left( {I^{2} + \left( {2a} \right)^{2}} \right)}.}}}} & (1) \end{matrix}$ Accordingly, the equations of motion for the element in the nth unit cell of the monoblock lattice are written as:

$\begin{matrix} {\mspace{79mu}{{m\frac{d^{\; 2}u_{n}}{{dt}^{2}}} = {K_{0}\left( {u_{n + 1} - {2u_{n}} + u_{n - 1}} \right)}}} & (2) \\ {\mspace{79mu}{{m\frac{d^{\; 2}v_{n}}{{dt}^{2}}} = {{K_{1}\left( {v_{n + 1} - {2v_{n}} + v_{n - 1}} \right)} + {\frac{{hK}_{1}}{2}\left( {\varphi_{n + 1} - \varphi_{n - 1}} \right)}}}} & (3) \\ {{I\frac{d^{\; 2}\varphi_{n}}{{dt}^{2}}} = {{K_{2}\left( {\varphi_{n + 1} - {2\varphi_{n}} + \varphi_{n - 1}} \right)} + {\frac{{hK}_{1}}{2}\left( {v_{n - 1} - v_{n + 1}} \right)} - {\frac{h^{2}K_{1}}{4}\left( {\varphi_{n + 1} + {2\varphi_{n}} + \varphi_{n - 1}} \right)}}} & (4) \end{matrix}$

The band structure of such a lattice would contain three bands corresponding to the three degrees of freedom (u, v, φ). However, for the sake of simplicity, this model is presently restricted to the propagation of rotational waves by allowing only rotation of the blocks about their center of mass and by constraining shear and longitudinal displacements in the lattice. The equation of motion associated with the rotational degrees of freedom (equation (4)) then takes the simpler general form:

$\begin{matrix} {{{I\frac{\partial^{\; 2}\varphi_{n}}{\partial t^{2}}} = {{K_{1}^{\prime}\left( {\varphi_{n + 1} - {2\varphi_{n}} + \varphi_{n - 1}} \right)} - {K_{2}^{\prime}\varphi_{n}}}}{where}{{K_{1}^{\prime} = {K_{2} - \frac{h^{2}K_{1}}{4}}},{{{and}\mspace{14mu} K_{2}^{\prime}} = {h^{2}{K_{1}.}}}}} & (5) \end{matrix}$

Dividing the equation by I yields the following rotational wave equation:

$\begin{matrix} {{{\frac{\partial^{\; 2}\varphi_{n}}{\partial t^{2}} - {\beta^{2}\left( {\varphi_{n + 1} - {2\varphi_{n}} + \varphi_{n - 1}} \right)} + {\alpha^{2}\varphi_{n}}} = 0}{with}{\beta^{2} = {{\frac{K_{1}^{\prime}}{I}\mspace{14mu}{and}\mspace{14mu}\alpha^{2}} = {\frac{K_{2}^{\prime}}{I}.}}}} & (6) \end{matrix}$

Equation (6) takes the form of the discrete Klein-Gordon equation. Equation (6) involves the second derivatives with respect to continuous time and the discrete second derivative with respect to position of the angular degree of freedom. Subsequently, following the approach of Dirac in linearizing the relativistic Klein-Gordon equation, a wave equation in terms of first order spatial and temporal derivatives can be derived by rewriting the Laplacian: Δφ_(n)=φ_(n+1)−2φ_(n)+φ_(n−1), in a “square root” form: Δφ_(n)=D(Dφ_(n)). This can be done exactly by introducing the following first order differential operator:

$\begin{matrix} {D = {{{e_{1}\Delta^{+}} + {e_{2}\Delta^{-}}} = {{\begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix}\Delta^{+}} + {\begin{pmatrix} 0 & 0 \\ 1 & 0 \end{pmatrix}\Delta^{-}}}}} & (7) \end{matrix}$ In equation (7), Δ⁺φ_(n)=φ_(n+1)−φ_(n) and Δ⁻φ_(n)=φ_(n)−φ_(n−1) are the forward and backward finite differences acting now on a two-vector. The 2×2 matrices e₁ and e₂ satisfy the conditions e₁e₂=e₂e₂=0 and e₁e₂+e₂e₁=I with I representing the 2×2 identity matrix. This formalism permits the exact and formal definition of the “square root” of the discrete Laplacian.

The Dirac-like equation for rotational waves corresponding to case I, then takes the form:

$\begin{matrix} {{\left\lbrack {{{\sigma_{x} \otimes I}\frac{\partial}{\partial t}} + {{i\;{{\beta\sigma}_{y} \otimes \left\{ {{e_{1}\Delta^{+}} + {e_{2}\Delta^{-}}} \right\}}} \pm {i\;\alpha\;{I \otimes I}}}} \right\rbrack\psi} = 0} & (8) \end{matrix}$

where σ_(x) and σ_(y) are the 2×2 Pauli matrices

${\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}\mspace{14mu}{and}\mspace{14mu}\begin{pmatrix} 0 & {- i} \\ i & 0 \end{pmatrix}},$ respectively. The parameter α plays the role of mass in the relativistic Dirac equation. Applying the outer product ⊗ leads to 4×4 matrices and ψ is a four-vector:

$\begin{pmatrix} \psi_{1n} \\ \psi_{2n} \\ \psi_{3n} \\ \psi_{4n} \end{pmatrix}.$

This four-component representation is the consequence of the discrete nature of the Laplacian. In contrast, with a continuous Laplacian, there is no distinction between forward and backward derivatives and only a two-component representation is needed. This is the case in the long-wavelength limit. In this limit, having a two-component spinor indicates that there is a coupling between waves propagating in opposite directions (positive or negative) along the chain of blocks.

When considering the short-wavelength four-component spinor solution, the first two components represent propagation of waves in the positive direction and the next two components propagation in the negative direction. The two components for the positive direction and the two components from the negative directions reflect a lifting of degeneracy due to asymmetry of the forward and backward finite different in the discrete Dirac equation. The solutions of equations (8) are automatically solutions of equation (6), but the converse is not true.

As will be seen later, the directions of propagation of the wave are expressed separately in the Dirac wave function. The ± in equation (8) corresponds to choices of the sign of the parameter α (i.e. choice of positive or negative “mass” in Dirac's formalism). First seeking solutions of equation (8) with the negative value, equation (8) becomes:

$\begin{matrix} {{\left\lbrack {{C\frac{\partial}{\partial t}} + {\beta\left\{ {{A\;\Delta^{+}} + {B\;\Delta^{-}}} \right\}} - {i\;\alpha\; I}} \right\rbrack\psi} = 0} & (9) \end{matrix}$ C, A, B, and I are the 4 matrices:

$\begin{pmatrix} 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \end{pmatrix},\begin{pmatrix} 0 & 0 & 0 & 1 \\ 0 & 0 & 0 & 0 \\ 0 & {- 1} & 0 & 0 \\ 0 & 0 & 0 & 0 \end{pmatrix},\begin{pmatrix} 0 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 \\ {- 1} & 0 & 0 & 0 \end{pmatrix},\begin{pmatrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{pmatrix}$

It is easily verifiable that CC=1, AA=BB=0, AB+BA=−1, and C (AΔ⁺+BΔ⁻)+(AΔ⁺+BΔ⁻) C=0, which are the conditions necessary to recover the wave equation (6) by applying the operator in equation (9) twice (with appropriate ±sign). Equation (9) is the basis for a discussion of rotational waves with non-trivial topologies in phononic structures.

Seeking solutions in the form of plane ψ_(jn)=a_(j)e^(−iωt)e^(iknh) with j=1, 2, 3, 4. ω and k are the angular frequency and wavenumber, respectively. Recalling that “h” is the spacing distance between blocks, equation (9) yields the system of equations:

$\begin{matrix} \left\{ \begin{matrix} {{{{- i}\;\alpha\; a_{1}} - {i\;\omega\; a_{3}} + {{\beta\left( {e^{ikh} - 1} \right)}a_{4}}} = 0} \\ {{{{- i}\;\alpha\; a_{2}} + {{\beta\left( {1 - e^{- {ikh}}} \right)}a_{3}} - {i\;\omega\; a_{4}}} = 0} \\ {{{{- i}\;\omega\; a_{1}} - {{\beta\left( {e^{ikh} - 1} \right)}a_{2}} - {i\;\alpha\; a_{3}}} = 0} \\ {{{{- {\beta\left( {1 - e^{- {ikh}}} \right)}}a_{1}} - {i\;\omega\; a_{2}} - {i\;\alpha\; a_{4}}} = 0} \end{matrix} \right. & (10) \end{matrix}$ This system of equations admits two doubly degenerate eigenvalues:

$\begin{matrix} {\omega = {{\pm \sqrt{\alpha^{2} - {{\beta^{2}\left( {e^{ikh} - 1} \right)}\left( {1 - e^{- {ikh}}} \right)}}} = {\pm \sqrt{\alpha^{2} + {4\beta^{2}\sin^{2}\frac{kh}{2}}}}}} & (11) \end{matrix}$

The negative frequency can be interpreted physically as follows. Since the angular field ought to be a real-valued quantity, it can be written as the sum of a complex term and its complex conjugate. The negative frequency is associated with the complex conjugate term. The rotational mode at the origin, k=0, has a finite frequency. Note that equation (11) gives two branches with positive and negative frequencies that do not intersect at the origin unless α=0.

The dispersion relations are periodic and defined in the first Brillouin zone:

$k \in {\left\lbrack {{- \frac{\pi}{h}},\frac{\pi}{h}} \right\rbrack.}$ Choosing the positive or negative branches of the dispersion relations, the four eigenvectors can be determined:

$\begin{matrix} {\psi_{1}^{\pm} = {{\begin{pmatrix} a_{1} \\ a_{2} \\ a_{3} \\ a_{4} \end{pmatrix}e^{{- i}\;\omega\; t}e^{iknh}} = {{a_{0}\begin{pmatrix} {{ie}^{{+ i}\frac{kh}{4}}\left( {\alpha - {\beta\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}} \right)} \\ {{\mp e^{{- i}\frac{kh}{4}}}\sqrt{\alpha^{2} - {\beta^{2}\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}^{2}}} \\ {{\mp {ie}^{{+ i}\frac{kh}{4}}}\sqrt{\alpha^{2} - {\beta^{2}\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}^{2}}} \\ {e^{{- i}\frac{kh}{4}}\left( {\alpha - {\beta\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}} \right)} \end{pmatrix}}e^{{- i}\;\omega\; t}e^{iknh}}}} & \left( {12a} \right) \\ {\psi_{2}^{\pm} = {{\begin{pmatrix} a_{1} \\ a_{2} \\ a_{3} \\ a_{4} \end{pmatrix}e^{{- i}\;\omega\; t}e^{iknh}} = {{a_{0}\begin{pmatrix} {- {{ie}^{{+ i}\frac{kh}{4}}\left( {\alpha + {\beta\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}} \right)}} \\ {{\mp e^{{- i}\frac{kh}{4}}}\sqrt{\alpha^{2} - {\beta^{2}\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}^{2}}} \\ {{\pm {ie}^{{+ i}\frac{kh}{4}}}\sqrt{\alpha^{2} - {\beta^{2}\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}^{2}}} \\ {e^{{- i}\frac{kh}{4}}\left( {\alpha + {\beta\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}} \right)} \end{pmatrix}}e^{{- i}\;\omega\; t}e^{iknh}}}} & \left( {12b} \right) \end{matrix}$

The upper signs in equations (12a) and (12b) correspond to the positive branch of the band structure (ω>0) and the lower signs to the negative branch (ω>0).

Similarly, the following equation can be solved:

$\begin{matrix} {{\left\lbrack {{C\frac{\partial}{\partial t}} + {\beta\left\{ {{A\;\Delta^{+}} + {B\;\Delta^{-}}} \right\}} + {i\;\alpha\; I}} \right\rbrack\overset{\_}{\psi}} = 0} & (13) \end{matrix}$

The operator in equation (13) is the complex conjugate of the operator in equation (9), so, seeking solutions in the form of plane waves, ψ _(jn)=σ _(j)e^(jωt)e^(−iknh) with j=1, 2, 3, 4 yields the four eigenvectors:

$\begin{matrix} {{\overset{\_}{\psi}}_{1}^{\pm} = {{\begin{pmatrix} {\overset{\_}{a}}_{1} \\ {\overset{\_}{a}}_{2} \\ {\overset{\_}{a}}_{3} \\ {\overset{\_}{a}}_{4} \end{pmatrix}e^{i\;\omega\; t}e^{- {iknh}}} = {{a_{0}\begin{pmatrix} {{ie}^{{- i}\frac{kh}{4}}\left( {\alpha - {\beta\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}} \right)} \\ {{\mp e^{{+ i}\frac{kh}{4}}}\sqrt{\alpha^{2} - {\beta^{2}\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}^{2}}} \\ {{\mp {ie}^{{- i}\frac{kh}{4}}}\sqrt{\alpha^{2} - {\beta^{2}\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}^{2}}} \\ {e^{{+ i}\frac{kh}{4}}\left( {\alpha - {\beta\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}} \right)} \end{pmatrix}}e^{i\;\omega\; t}e^{- {iknh}}}}} & \left( {14a} \right) \\ {{\overset{\_}{\psi}}_{2}^{\pm} = {{\begin{pmatrix} {\overset{\_}{a}}_{1} \\ {\overset{\_}{a}}_{2} \\ {\overset{\_}{a}}_{3} \\ {\overset{\_}{a}}_{4} \end{pmatrix}e^{i\;\omega\; t}e^{- {iknh}}} = {{a_{0}\begin{pmatrix} {- {{ie}^{{- i}\frac{kh}{4}}\left( {\alpha + {\beta\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}} \right)}} \\ {{\mp e^{{+ i}\frac{kh}{4}}}\sqrt{\alpha^{2} - {\beta^{2}\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}^{2}}} \\ {{\pm {ie}^{{- i}\frac{kh}{4}}}\sqrt{\alpha^{2} - {\beta^{2}\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}^{2}}} \\ {e^{{+ i}\frac{kh}{4}}\left( {\alpha + {\beta\left( {e^{i\frac{kh}{2}} - e^{{- i}\frac{kh}{2}}} \right)}} \right)} \end{pmatrix}}e^{i\;\omega\; t}e^{- {iknh}}}}} & \left( {14b} \right) \end{matrix}$ Non-Conventional Topology of Rotational Waves

Turning to an investigation of the topology in k space of the solutions given by equations (12a), (12b) and (14a), (14b), it is easy to show the following relations:

$\begin{matrix} {{\psi_{1,2}^{\pm}\left( {k + \frac{2\pi}{h}} \right)} = {e^{{- i}\frac{\pi}{2}}{{\overset{\_}{\psi}}_{1,2}^{* \pm}(k)}}} & \left( {15a} \right) \\ {{\psi_{1,2}^{\pm}\left( {k + \frac{4\pi}{h}} \right)} = {e^{{- i}\;\pi}{\psi_{1,2}^{\pm}(k)}}} & \left( {15b} \right) \\ {{\psi_{1,2}^{\pm}\left( {k + \frac{6\pi}{h}} \right)} = {e^{{- i}\frac{3\pi}{2}}{{\overset{\_}{\psi}}_{1,2}^{* \pm}(k)}}} & \left( {15c} \right) \\ {{\psi_{1,2}^{\pm}\left( {k + \frac{8\pi}{h}} \right)} = {e^{{- i}\; 2\;\pi}{\psi_{1,2}^{\pm}(k)}}} & \left( {15d} \right) \end{matrix}$

In equations (15a)(15d), the * symbol represents the complex conjugate. ψ and ψ are orthogonal and satisfy: ψψ=0. The relations given by equation (15) indicate that the wave function has a periodicity of

$\frac{8\pi}{h}..$ The topology of the wave function is more easily illustrated by considering the parallel transport of two orthogonal vector fields along a closed loop corresponding to

${k \in \left\lbrack {0,\frac{8\pi}{h}} \right\rbrack}..$ The first vector field represents the wave function ψ and the second vector field the complex conjugate of ψ. This non-conventional torsional topology is illustrated in FIG. 2 in the form of a torus 200 with a square cross section 205 that exhibits four 90° twists.

Starting at the top of the figure, the red arrow and green arrow represent the two orthogonal wave functions ψ and ψ respectively. The red arrow points in a direction perpendicular to the plane of the figure and the green arrow points in a direction within the plane of the figure. As the vector fields are transported along the first twist, the green arrow is now pointing in the plane of the figure and the green arrow becomes perpendicular to that plane, that is ψ turns into ψ.

The red arrow transported through the first twist now makes a right angle with the red arrow at the starting point. This means that the wave function ψ has accumulated a phase of

$\frac{- \pi}{2}.$ The process is repeated through the second twist. The red arrow transported through the second twist is now parallel to the red arrow at the origin but is pointing in the opposite direction and satisfies the condition:

${\psi\left( {k + \frac{4\;\pi}{h}} \right)} = {- {{\overset{\_}{\psi}(k)}.}}$ This means that the wave function ψ has accumulated a phase of −π. The vector fields can be transported through the third and fourth twists (not illustrated in the figure) to close the loop in k space and reach the starting point. At that point:

${\psi\left( {k + \frac{8\pi}{h}} \right)} = {{\psi(k)}.}$

The torsional topology of the wave function results from its spinor characteristics that endow it with fermion-like character. This character is analyzed in further detail below. To simplify the mathematics, only the long-wavelength limit is considered.

Analysis of Fermion-Like Behavior in the Long-Wavelength Limit

Dirac Equation in the Long-Wavelength Limit

In the long-wavelength limit, k→0, equation (6) reduces to

$\begin{matrix} {{\frac{\partial^{2}\varphi}{\partial t^{2}} - {\beta^{2}\frac{\partial^{2}\varphi}{\partial x^{2}}} + {\alpha^{2}\varphi}} = 0} & (16) \end{matrix}$

Note that equation (16) is also that of a system composed of a one-dimensional harmonic chain 301 composed of masses 302 (m) connected to their nearest neighbors by springs 303 with spring constant k₁ and connected to a rigid substrate by springs 304 which constant is k₂ (see FIG. 3). The relation between the equation of this spring-mass system 300 and the rotational wave equation (16) is obtained by setting α²κ₂/m and β²=κ₁/m. In the long-wavelength limit, this is a system composed of a string embedded in an elastic medium. The isomorphism between the block-spring system and the mass-spring model will be used in subsequent calculations.

Equation (16) admits plane wave solutions with eigenvalues given by ω=±√{square root over (α²+β²(kh)²)}. At k=0, the eigenmode with frequency ω=α correspond to point B in FIG. 2. Note that to obtain the dispersion relation, equation (16) was solved by assuming plane wave solutions. These solutions correspond to traveling waves. However, the dispersion relation indicates that, for instance, at k=0 the group velocity is zero, that is the eigenvector solution to (16) is a standing wave. Since a standing wave is the superposition of two waves traveling in opposite directions, equation (16) is unable to reveal the standing wave character of its solutions. Furthermore, the amplitude of the two oppositely traveling waves forming a standing wave have to be the same. Equation (16) does not impose any constraint on the amplitude of its plane wave solutions. For k≠0, the group velocity approaches the phase velocity only when k→∞, the solutions of equation (16) have mixed standing and traveling wave character. These solutions ought to be the superposition of waves traveling in opposite directions, whose amplitudes are not independent of each other, but are wavenumber dependent. It is these very observations that make the factorization of equation (16) into direction-specific equations necessary.

Equation (16) is isomorphic to the continuous relativistic Klein-Gordon equation. Here, again, following the approach of Dirac in linearizing the relativistic Klein Gordon equation, the goal is to derive a wave equation in terms of first order spatial and temporal derivatives:

$\begin{matrix} {{\left\lbrack {{\sigma_{X}\frac{\partial}{\partial t}} + {i\;\beta\;\sigma_{Y}\frac{\partial}{\partial_{X}}} - {i\;\alpha\; I}} \right\rbrack\psi} = 0} & \left( {17a} \right) \\ {{\left\lbrack {{\sigma_{X}\frac{\partial}{\partial t}} + {i\;\beta\;\sigma_{Y}\frac{\partial}{\partial_{X}}} - {i\;\alpha\; I}} \right\rbrack\overset{\_}{\psi}} = 0} & \left( {17b} \right) \end{matrix}$ where σ_(X) and σ_(Y) are the 2×2 Pauli matrices and I is the 2×2 identity matrix.

Note that these equations are not self-dual. Taking the complex conjugate of equation (17a) results in equation (17b), indeed the first two terms are real and only the last term changes sign. In particular, this results from the negative sign of the second term in the Klein-Gordon equation (16), which involves the multiplicative imaginary number “I” in the second term of the Dirac-like equations (17a), (17b).

Then ψ=ψ. S_(D) while ψ is a solution of equation (17a), its complex conjugate is not a solution to (17a). ψ is solution to equation (17b). In the language of Quantum Field Theory, ψ and ψ represent two different physical entitles, namely “particles” and “antiparticles.” 3.2 Lagrangian Formalism

The Dirac-like equation (17a) derived in the previous section can be obtained from the Lagrangian:

$\begin{matrix} {{L\left( {\overset{\_}{\psi},\frac{\partial\overset{\_}{\psi}}{\partial t},\frac{\partial\overset{\_}{\psi}}{\partial x}} \right)} = {{\frac{\partial\overset{\_}{\psi}}{\partial t}\sigma_{x}\psi} + {i\;\beta\frac{\partial\overset{\_}{\psi}}{\partial X}\sigma_{y}\psi} + {i\;\alpha\overset{\_}{\psi}\; I\;\psi}}} & \left( {18a} \right) \end{matrix}$ In equation (18a), ψ is the Hermitian conjugate of ψ. The Lagrangian (18a) is Hermitian and its conjugate is given by:

$\begin{matrix} {{\overset{\_}{L}\left( {\psi,\frac{\partial\psi}{\partial t},\frac{\partial\psi}{\partial x}} \right)} = {{\overset{\_}{\psi}\sigma_{x}\frac{\partial\psi}{\partial t}} + {i\;\beta\overset{\_}{\psi}\sigma_{y}\frac{\partial\overset{\_}{\psi}}{\partial X}} - {i\;{\alpha\psi}\; I\;\overset{\_}{\psi}}}} & \left( {18b} \right) \end{matrix}$ The Lagrange's equation motion for the field variables ψ is given by:

$\begin{matrix} {{{\frac{\partial}{\partial t}\left( \frac{\partial L}{\partial\left( {\partial_{t}\left( \overset{\_}{\psi} \right)} \right.} \right)} + {\frac{\partial}{\partial x}\left( \frac{\partial L}{\partial\left( {\partial_{x}\left( \overset{\_}{\psi} \right)} \right.} \right)} - \frac{\partial L}{\partial_{x}\left( \overset{\_}{\psi} \right)}} = 0} & (19) \end{matrix}$ An equation of motion for the field ψ can also be obtain from equation (19) by replacing L by L and ψ by ψ. Using this approach, it is straightforward to recover equations (17a) and (17b). Solutions of the Long-Wavelength Dirac Equation

Having established the Dirac equation and its Hermitian conjugate (17a), (17b), it is expected that there exist solutions in the form of plane waves, and since these are matrix equations, there exists a solution that contains a spinor part (i.e. a two by one part matrix). The solutions may therefore be written to be presented in the form: ψ_(k)=ψ(k,ω _(k))=c ₀ξ_(k))e ^((±)iω) ^(k) ^(t) e ^((±)ikx) and ψ _(k)=ψ(k,ω _(k))=c ₀ ξ _(k))e ^((±)iω) ^(k) ^(t) e ^((±)ikx) where ξ_(k) and ξ _(k) are two by one spinors. Inserting the various forms for these solutions in equations (17a), (17b) lead to the same eigenvalues that were obtained with the Klein-Gordon equation, namely: ω=±√{square root over (α²+β²(kh)²)}. The spinor part the solutions for the different plane waves is summarized in Table 1.

TABLE 1 Two by one spinor solutions of equations (17a) and (17b) for the different plane waveforms. e^(+ikx)e^(+iω) ^(k) ^(t) e^(−ikx)e^(+iω) ^(k) ^(t) e^(+ikx)e^(+iω) ^(k) ^(t) e^(−ikx)e^(+iω) ^(k) ^(t) ξ_(k) $\quad\begin{pmatrix} \sqrt{\omega + {\beta kh}} \\ \sqrt{{\omega - {\beta kh}}\mspace{11mu}} \end{pmatrix}$ $\quad\begin{pmatrix} \sqrt{{\omega - {\beta kh}}\;} \\ \sqrt{{\omega + {\beta kh}}\;} \end{pmatrix}$ $\quad\begin{pmatrix} {- \sqrt{\omega - {\beta kh}}} \\ \sqrt{{\omega + {\beta kh}}\;} \end{pmatrix}$ $\quad\begin{pmatrix} {- \sqrt{\omega + {\beta kh}}} \\ \sqrt{{\omega - {\beta kh}}\;} \end{pmatrix}$ ξ _(k) $\quad\begin{pmatrix} \sqrt{\omega - {\beta\;{kh}}} \\ {- \sqrt{{\omega + {\beta\;{kh}}}\;}} \end{pmatrix}$ $\quad\begin{pmatrix} \sqrt{\omega + {\beta\;{kh}}} \\ {- \sqrt{{\omega - {\beta\;{kh}}}\;}} \end{pmatrix}$ $\quad\begin{pmatrix} \sqrt{\omega + {\beta\;{kh}}} \\ \sqrt{{\omega - {\beta\;{kh}}}\;} \end{pmatrix}$ $\quad\begin{pmatrix} \sqrt{\omega - {\beta\;{kh}}} \\ \sqrt{{\omega + {\beta\;{kh}}}\;} \end{pmatrix}$

The spinors in Table 1 are given to within any complex constant c₀.

Note that the wave functions ψ ψ will eventually be promoted to the rank of creation and annihilation operators. Additionally, if α=0, equation (17a) becomes

${\left\lbrack {{\sigma_{x}\frac{\partial}{\partial t}} + {i\;{\beta\sigma}_{y}\frac{\partial}{\partial x}}} \right\rbrack\psi} = 0.$ Using a plane wave solution with

${{\xi_{k}\left( {k,\omega_{k}} \right)} = \begin{pmatrix} \alpha_{1} \\ \alpha_{2} \end{pmatrix}},$ this equation reduces to the system:

$\left\{ {\begin{matrix} {{\left( {\omega - {\beta\; k\; h}} \right)\alpha_{1}} = 0} \\ {{\left( {\omega + {\beta\; k\; h}} \right)\alpha_{2}} = 0} \end{matrix}.} \right.$

Two solutions can be obtained for the angular velocity of the plane wave ω=±βkh. These correspond to plane waves propagating in the positive and negative directions. In this case, the components of the two-spinor, a₁ and a₂ are now independent of each other and independent of the wavenumber. The amplitude of the plane wave propagating in the positive direction is independent of that of the wave propagating in the opposite direction. When α≠0, from Table 1, the components of the two-spinor are not independent of each other. This indicates that the directions of propagation are not independent of each other anymore. It is the parameter α that couples those directions.

Energy and Anticommutation

The Hamiltonian density associated with the Lagrangian, L. It is given by the Legendre transformation:

$\begin{matrix} {\overset{\_}{H} = {{\overset{\_}{\Pi}\frac{\partial\psi}{\partial t}} - \overset{\_}{L}}} & (20) \end{matrix}$ Where the momentum conjugate:

$\overset{\_}{\Pi} = {\frac{\partial\overset{\_}{L}}{\partial\left( {{\partial\psi}/{\partial t}} \right)} = {\overset{\_}{\psi}{\sigma_{x}.}}}$ Inserting this expression in Equation (20) and using equation (18b), yields the Hamiltonian density:

$\begin{matrix} {\overset{\_}{H} = {\overset{\_}{\psi}\;\sigma_{x}\frac{\partial\psi}{\partial t}}} & (21) \end{matrix}$

As mentioned previously, the quantity ψ is now promoted to an operator (using the same symbol ψ to represent the operator) and a representation is used whereby the field ψ is expanded into plane waves and their complex conjugate:

$\begin{matrix} {{\psi(x)} = {\sum\limits_{k}\;{\frac{1}{\sqrt{2\;\omega}}\left\lbrack {{\alpha_{k}\xi_{k}e^{i\; k\; x}e^{{- i}\;\omega\; t}} + {\alpha_{k}^{*}\xi_{k}^{*}e^{{- i}\; k\; x}e^{i\;\omega\; t}}} \right\rbrack}}} & (22) \end{matrix}$

Similarly, the Hermitian conjugate of the ψ operator can be written:

$\begin{matrix} {{\overset{\_}{\psi}(x)} = {\sum\limits_{k}\;{\frac{1}{\sqrt{2\;\omega}}\left\lbrack {{{\overset{\_}{\alpha}}_{k}{\overset{\_}{\xi}}_{k}e^{{- i}\; k\; x}e^{i\;\omega\; t}} + {{\overset{\_}{\alpha}}_{k}^{*}{\overset{\_}{\xi}}_{k}^{*}e^{i\; k\; x}e^{{- i}\;\omega\; t}}} \right\rbrack}}} & (23) \end{matrix}$ The quantities ā_(k), ā_(k) ^(k), a_(k), and a_(k) ^(k) are creation and annihilation operators, respectively. Inserting equations (22) and (23) into equation (21) and integrating over all space results in:

$\begin{matrix} {E = {{\int{{dx}\overset{\_}{H}}} = {{\sum\limits_{k}{\sum\limits_{k^{\prime}}{\frac{{\overset{\_}{a}}_{k^{\prime}}}{\sqrt{2\omega^{\prime}}}\left( {{- i}\;\omega} \right){\overset{\_}{\xi}}_{k^{\prime}}\sigma_{x}\xi_{k}\frac{a_{k}}{\sqrt{2\omega}}e^{{- {i{({\omega - \omega^{\prime}})}}}t}{\int{dxe}^{{- {i{({k^{\prime} - k})}}}x}}}}} + {\sum\limits_{k}{\sum\limits_{k^{\prime}}{\frac{{\overset{\_}{a}}_{k^{\prime}}}{\sqrt{2\omega^{\prime}}}\left( {i\;\omega} \right){\overset{\_}{\xi}}_{k^{\prime}}\sigma_{x}\xi_{k}^{*}\frac{a_{k}^{*}}{\sqrt{2\omega}}e^{{i{({\omega + \omega^{\prime}})}}t}{\int{dxe}^{{- {i{({k^{\prime} + k})}}}x}}}}} + {\sum\limits_{k}{\sum\limits_{k^{\prime}}{\frac{{\overset{\_}{a}}_{k^{\prime}}^{*}}{\sqrt{2\omega^{\prime}}}\left( {{- i}\;\omega} \right){\overset{\_}{\xi}}_{k^{\prime}}^{*}\sigma_{x}\xi_{k}\frac{a_{k}}{\sqrt{2\omega}}e^{{- {i{({\omega + \omega^{\prime}})}}}t}{\int{dxe}^{{i{({k^{\prime} + k})}}x}}}}} + {\sum\limits_{k}{\sum\limits_{k^{\prime}}{\frac{{\overset{\_}{a}}_{k^{\prime}}^{*}}{\sqrt{2\omega^{\prime}}}\left( {i\;\omega} \right){\overset{\_}{\xi}}_{k^{\prime}}^{*}\sigma_{x}\xi_{k}^{*}\frac{a_{k}^{*}}{\sqrt{2\omega}}e^{{- {i{({\omega^{\prime} - \omega})}}}t}{\int{dxe}^{{- {i{({k - k^{\prime}})}}}x}}}}}}}} & (24) \end{matrix}$

The spatial integrals in equation (24) can be replaced by their respective delta functions, δ_(k′,k), δ_(k′,−k). Using Table 1, it can be shown that ξ _(k)σ_(k)ξ_(k)=ξ _(k)*σ_(k)ξ_(k)*=2ωc₀ ² and it can be shown that ξ _(−k)σ_(k)ξ_(k)*=ξ _(−k)σ_(k)ξ_(k)=0. Taking the constant c₀ ²=−i to obtain an energy in the form of a real number, the energy of the system becomes:

$\begin{matrix} {E = {\sum\limits_{k}{\omega\left( {{{\overset{\_}{a}}_{k}a_{k}} - {{\overset{\_}{a}}_{k}^{*}a_{k}^{*}}} \right)}}} & (25) \end{matrix}$

To avoid the negative contribution to the energy in equation (25) and following the methods of quantum field theory, anticommutation rules need to be imposed on the creation and annihilation operators. That is: [a _(k),σ_(κ)]=ā _(k) a _(κ) +a _(κ) ā _(k)=δ_(k,κ) With these anticommutation rules, the energy takes the form

$\begin{matrix} {{E - E_{0}} = {\sum\limits_{k}{\omega\left( {{{\overset{\_}{a}}_{k}a_{k}} + {a_{k}^{*}{\overset{\_}{a}}_{k}^{*}}} \right)}}} & (27) \end{matrix}$

In equation (27), the terms that arise from the delta function of equation (26) have been lumped into a negative zero-point energy E₀=−Σ_(k)ω(k).

The anticommutation rules given by equation (26) also imply that ā_(k)kā_(k)=0. This is Pauli's Exclusion Principle, which states that more than one particle per state cannot be created. The anticommutation rules indicate that the rotational modes have fermion character. A physical interpretation of these anticommutation rules suggests that the fields given by (22) and (23) are composed of forward and backward traveling waves whose relative amplitudes are determined solely by the spinor part of the wave function ξ_(k)e^(ikx)e^(−iωt). Attempting to add an additional traveling component to one of the quasi-standing solutions via the operator a_(k)*ā_(k)* would create a new state that could not be a solution of the Dirac equations (17a), (17b). Such a state is not allowed. This restriction provides a physical interpretation of Pauli's Exclusion Principle.

Number Operators

It was shown above that the energy is an invariant. This invariant is the consequence of the time reversal symmetry of equations (17a) and (17b). There exist other symmetries. In particular the Dirac-like equations are invariant upon the transformation involving a phase shift. Consequently, it can be shown that the number operator: N=∫dxψσ _(x)ψ is also an invariant. Starting from the field operator given by equation (22), k can be replaced by −k in the second term of the summation. This leads to:

$\begin{matrix} {{\psi(x)} = {{\sum\limits_{k}{\frac{1}{\sqrt{2\omega}}\left\lbrack {{a_{k}\xi_{k}e^{ikx}e^{{- i}\;\omega\; t}} + {a_{- k}^{*}\xi_{- k}^{*}e^{ikx}e^{{- i}\;\omega\; t}}} \right\rbrack}} = {\sum\limits_{k}{\frac{1}{\sqrt{2\omega}}\left\lbrack {{a_{k}\xi_{k}e^{ikx}e^{{- i}\;\omega\; t}} + {b_{k}\eta_{k}e^{ikx}e^{i\;\omega\; t}}} \right\rbrack}}}} & (29) \end{matrix}$

Similarly, the following field can be written:

$\begin{matrix} {{\overset{\_}{\psi}(x)} = {{\sum\limits_{k}{\frac{1}{\sqrt{2\omega}}\left\lbrack {{{\overset{\_}{a}}_{k}{\overset{\_}{\xi}}_{k}e^{- {ikx}}e^{i\;\omega\; t}} + {{\overset{\_}{a}}_{- k}^{*}{\overset{\_}{\xi}}_{- k}^{*}e^{- {ikx}}e^{{- i}\;\omega\; t}}} \right\rbrack}} = {\sum\limits_{k}{\frac{1}{\sqrt{2\omega}}\left\lbrack {{{\overset{\_}{a}}_{k}{\overset{\_}{\xi}}_{k}e^{- {ikx}}e^{i\;\omega\; t}} + {{\overset{\_}{b}}_{k}{\overset{\_}{\eta}}_{k}e^{- {ikx}}e^{{- i}\;\omega\; t}}} \right\rbrack}}}} & (30) \end{matrix}$

Inserting equations (29) and (30) into (28), and following the same approach as before, yields:

$\begin{matrix} {N = {\sum\limits_{k}{\frac{1}{2\omega}\left( {{{\overset{\_}{a}}_{k}a_{k}{\overset{\_}{\xi}}_{k}\sigma_{x}\xi_{k}} + {{\overset{\_}{b}}_{k}b_{k}{\overset{\_}{\eta}}_{k}\sigma_{x}\eta_{k}}} \right)}}} & (31) \end{matrix}$

Using Table 1 for the spinors and their Hermitian conjugate yields:

$\begin{matrix} {N = {\sum\limits_{k}\left( {{{\overset{\_}{a}}_{k}a_{k}} + {{\overset{\_}{b}}_{k}b_{k}}} \right)}} & (32) \end{matrix}$ This is clearly an invariant. The operator ā_(k)a_(k) count the number of waves propagating in a given direction and b _(k)b_(k) correspond to the number of waves propagating in the opposite direction. Considering a field representation that includes only one type of plane wave and therefore, looking at the first term of equation (32), provides a definition for the operator:

$\begin{matrix} {N_{k} = {\frac{1}{2\omega}{\overset{\_}{a}}_{k}a_{k}{\overset{\_}{\xi}}_{k}\sigma_{x}\xi_{k}}} & (33) \end{matrix}$

Recalling that ξ _(k)σ_(k)ξ_(k)=2ω=ω+βkh+ω−βkh. the spinors can be renormalized with respect to the arbitrary complex constant c₀. The last equality suggests the introduction of direction switching operators

$S_{+} = {{\frac{1}{2}\left( {\sigma_{x} + {i\;\sigma_{y}}} \right)} = {{\begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix}\mspace{14mu}{and}\mspace{14mu} S_{-}} = {{\frac{1}{2}\left( {\sigma_{x} - {i\;\sigma_{y}}} \right)} = {\begin{pmatrix} 0 & 0 \\ 1 & 0 \end{pmatrix}.}}}}$ These operators anticommute: S₊S⁻+S⁻S₊=I. Using this later relation permits equation (33) to be reformulated in the form:

$\begin{matrix} {N_{k} = {{\frac{1}{2\omega}{\overset{\_}{a}}_{k}a_{k}{\overset{\_}{\xi}}_{k}I\;\sigma_{x}\xi_{k}} = {{\overset{\_}{a}}_{k}a_{k}{{\overset{\_}{\xi}}_{k}\left( {{\frac{1}{2\omega}S_{+}S_{-}\sigma_{x}} + {\frac{1}{2\omega}S_{-}S_{+}\sigma_{x}}} \right)}\xi_{k}}}} & (34) \end{matrix}$ The operator

${\frac{1}{2\omega}S} + S - \sigma_{x}$ defines the occupancy of one of the directions along the chain of mass and springs. Its eigenvalues are given by

$n_{k}^{+} = {{{{\overset{\_}{\xi}}_{k}\frac{1}{2\omega}S} + S - {\sigma_{x}\xi_{k}}} = {\frac{\omega + {\beta\;{kh}}}{2\omega}.}}$ The operator

${\frac{1}{2\omega}S} - S + \sigma_{x}$ defines the occupancy of the other opposite direction. Its eigenvalues are given by:

$n_{k}^{-} = {{{{\overset{\_}{\xi}}_{k}\frac{1}{2\omega}S} - S + {\sigma_{x}\xi_{k}}} = {\frac{\omega - {\beta\;{kh}}}{2\omega}.}}$ Considering the dispersion relation for the positive branch of the band structure, ω=+√{square root over (α²+β²(kh)²)}, the eigenvalues n_(k) ⁺ and n_(k) ⁻, can be rewritten in terms of the wave vector k only:

$\begin{matrix} {n_{k}^{\pm} = {\frac{1}{2} \pm \frac{\beta\;{{kh}/\alpha}}{2\sqrt{1 + \left( \frac{\beta\;{kh}}{\alpha} \right)^{2}}}}} & (35) \end{matrix}$ At the origin of the band structure, k=0, and n_(k) ⁺=n_(k) ⁻=½. The direction occupancy is the same for the two opposite directions; hence the state of the system is described by a standing wave. For k→∞, n_(k) ⁺=1, and n_(k) ⁻=0, the phononic structure supports a pure traveling wave with only one direction favored. For finite k, one has the superposition of traveling and standing waves. By spanning the dispersion relation, one spans the “direction” states. One can quantify the “amount” of traveling wave character of the wave function with the quantity: (n_(k) ⁺−n_(k) ⁻). When expressed as a function of frequency, this quantity takes the functional form:

$\left( {n_{k}^{+} - n_{k}^{-}} \right) = {\frac{\sqrt{\omega^{2} - \alpha^{2}}}{\omega}.}$ This clearly shows that the wave function character evolves from a pure standing wave at ω=α to a pure traveling wave as the frequency increases to infinity.

The fermion-like behavior of the rotational wave opens up opportunities in the control of the direction of propagation of elastic waves. For instance, it has been shown that the two-spinor solution represents coupling between waves propagating in opposite directions. The spinor part of the wave function can be projected on the orthonormal basis

$\begin{pmatrix} 1 \\ 0 \end{pmatrix}\mspace{14mu}{and}\mspace{14mu}\begin{pmatrix} 0 \\ 1 \end{pmatrix}$ representing the possible directions of propagation of the wave. This enables information to be encoded in the relative weight (phase) of the directions of propagation by controlling the wavenumber, k. This control can enable applications in the field of information processing, as discussed below. Applications to Information Processing Phase Control (ψ)-Bit

The properties of rotational waves and their fermion-like behavior can be used to encode and process information, relying upon the isomorphism between the block-spring model and a mass-spring system (see FIG. 3) when considering a device that is composed of a segment of mass-spring system isomorphic to the block-spring system supporting the rotational waves sandwiched between two semi-infinite harmonic chains (FIG. 4).

The number of masses in the central region of the sandwich is denoted N_(c). Assume that all the masses, m, have the same value, and that the spring constants of the semi-infinite regions 401 are the same as the mass-to-mass spring constant, K₁. The spring constant of the side springs 404 is K₂. Recall that these parameters map onto the parameters characterizing the rotational wave equation via the relations: α²=K₂/m and β²=K₁/m.

Of particular interest is the transmission of plane waves launched from the left semi-infinite chain through the fermion-like segment. For this, the transmission coefficient of the composite system 400 of FIG. 4 is calculated, employing the methods of the Interface Response Theory (IRT). The use of this method to calculate the transmission coefficient of infinite harmonic chains 401 perturbed with harmonic finite side chains 404 is known to those of ordinary skill in the art.

FIG. 5 depicts the transmission coefficient as a function of frequency for a segment composed of N_(c)=80 masses with side springs, with m=1 and K₁=1, i.e. β²=1. The side springs 404 have stiffness K₂=0.1, i.e. α=0.316.

The transmission coefficient is zero for frequencies up to the value ω=α=0.316. This zero transmission corresponds to the gap in the band structure of the rotational waves. Waves that are launched and propagate in the right semi-infinite harmonic chain cannot propagate through the central segment of the composite. Inside the passing band of the central segment, transmission occurs. The oscillatory nature of the transmission coefficient is due to the finite size and the discrete nature of the central segment. The frequency equal to 2 corresponds to the upper frequency of waves that can be supported by a discrete one-dimensional harmonic chain. The decrease of the transmission coefficient as the frequency of 2 is approached results from the curvature (zero group velocity) of the band of the one-dimensional harmonic chain. Therefore, a one-to-one correspondence between the long-wavelength limit result (dotted line) and the transmission through a discrete system cannot be established.

However, the general trend of an increase in transmission as the frequency increases from the bottom of the passing band of the central segment structure is representative of both the discrete and long-wavelength limit fermion-like behavior. The most important observation is that a one-to-one correspondence between a measurable scalar quantity can be established: the transmission coefficient and the spinor characteristics of the wave that can be supported by the central segment of the composite.

For instance, for an incident wave e^(ikx)e^(iωt) at ω=α (k=0), the transmission coefficient is zero and the spinor part of the wave function takes on the form:

${\alpha\begin{pmatrix} 1 \\ 1 \end{pmatrix}}.$ The components of the spinor vary monotonically as the frequency is increased (i.e. the wavenumber k is increased) until it takes a form approaching:

$\sqrt{2\beta\;{kh}}{\begin{pmatrix} 1 \\ 0 \end{pmatrix}.}$ A similar tuning of the spinor part of the wave function can be achieved by fixing the frequency of the incident waves and varying the stiffness κ₂ (i.e. α) of the side springs relative to the stiffness of the 1D harmonic chain. From a practical point of view, this could be achieved by using springs composed of an elastic material whose stiffness is dependent on the magnitude of an externally applied electromagnetic wave (e.g., materials exhibiting photo-elastic effects) or on an external magnetic field (e.g., a magneto-elastic medium).

Defining a basis set |0) and |1) whose matrix representation takes the form of

$\begin{pmatrix} 1 \\ 0 \end{pmatrix}\mspace{14mu}{and}\mspace{14mu}\begin{pmatrix} 0 \\ 1 \end{pmatrix}$ permits me creation of a superposition of state u(ω)|0)+v(ω)|1) at fixed α or u(α)|0)+v(α)|1) at fixed ω. Note that launching the incident wave from the right-hand side of the composite harmonic chain corresponds to an incident wave with ω<0. In that case, for instance for ω=α (k=0), the spinor takes the value

${\alpha\begin{pmatrix} {- 1} \\ 1 \end{pmatrix}}.$

This spinor is equivalent to the value

${\alpha\begin{pmatrix} 1 \\ {- 1} \end{pmatrix}}\quad$ to within a factor of −1. Therefore, for a given incident wave of frequency ω, the operation of tuning α can be visualized between zero and ω=α as that of a Hadamard transformation. As seen previously, the Hadamard transformation can then be physically achieved by modulating the strength of an external stimulus such as the magnitude of an external field. In its matrix representation, the Hadamard transformation takes the form

$\begin{pmatrix} 1 & 1 \\ 1 & {- 1} \end{pmatrix}\quad$ and transforms a wave launched from the left

$\begin{pmatrix} 1 \\ 0 \end{pmatrix}\quad$ into the standing wave, whose spinor is proportional to

$\begin{pmatrix} 1 \\ 1 \end{pmatrix}{\quad.}$ It also transforms a wave propagating from the right

$\begin{pmatrix} 0 \\ 1 \end{pmatrix}\quad$ into me standing wave whose spinor part is proportional to

$\begin{pmatrix} 1 \\ {- 1} \end{pmatrix}.$

The action of applying the Hadamard gate a second time is equivalent to detuning the parameter α between from ω=α to zero. This transforms the states

$\begin{pmatrix} 1 \\ 1 \end{pmatrix}\mspace{14mu}{and}\mspace{14mu}\begin{pmatrix} 1 \\ {- 1} \end{pmatrix}$ back into the traveling states

$\begin{pmatrix} 1 \\ 0 \end{pmatrix}\mspace{14mu}{and}\mspace{14mu}{\begin{pmatrix} 0 \\ 1 \end{pmatrix}.}$ This transformation can be detected easily through the change in transmission from zero to 1. The control of the relative magnitude and sign of the spinor components of the wave function enables the encoding of information in the phase of the superposition of states. A system like that of FIG. 4 is defined as a ψ(phase)-bit. The ψ-bit can be operated on through the system's physical parameter α via some external stimulus to transform the spinor state of the wave function. The Hadamard gate is an example of such an operation. The spinor state of the central segment of the ψ-bit can be determined by measuring a single scalar quantity, namely the transmission coefficient. The question that then arises is that of the encoding of information on multiple ψ-bits and the possibility of operating on them, addressed below. Operations on Multiple ψ-Bits

Consider ψ-bits, i=1, 2, each prepared in the state: |ξ_(i)

=u(α_(i))|0

+v(α_(i))|1

. The side springs in these two ψ-bits are again composed of material whose stiffness is dependent on some external stimulus. Two external stimuli are applied on each ψ-bits independently at this time. As seen in the previous section, this is achieved by launching a wave with frequency ω>α_(i) to the right in each left semi-infinite chain. The transmission coefficient is measurable for each bit and takes on the value T(α_(i)). The two-bit system is in the initial product state: |ξ₁)⊗|ξ₂). One can now operate on this product state by defining a function ƒ(x) such that the product state becomes: |ξ₁)⊗ƒ(ξ₁)⊕ξ₂). Because there is a one-to-one relation between T(α₁) and |ξ₁), this can be achieved by constructing a physical device that links the value of α₂ (i.e., the magnitude of the external stimulus acting on the second ψ-bit) to that of the transmission coefficient of the first ψ-bit. For instance, one can define α₂(T (α₁))=ƒ(T(α₁)). The new state of the second ψ-bit will therefore take the form |ξ₂)=u(α₂(T(α₁)))|0)+v(α₂(T(α₁)))|1). This procedure enables parallelism in the application of the operating function. This type of operation can be applied straightforwardly to a system of N ψ-bits, which could lead to parallel operations on the N-bit product state.

CONCLUSIONS

A 1D discrete block-spring model that supports rotational waves was investigated. This 1D discrete model was then analyzed within the Dirac's formalism to show that rotational waves are characterized by wave functions that include a spinor part and a spatio-temporal part. The spinor part is related to the coupling between the senses of wave propagation. The spatio-temporal part of the wave function retains its plane-wave character.

The topology of the wave function in k-space was also investigated to show that the wave function possesses a non-conventional torsional topology. The wave function is supported by a square cross section torus manifold exhibiting four 90° twists. This manifold has a periodicity of

$\frac{g_{\pi}}{h},$ where h is the periodicity of the 1D phononic structure.

The 1D system, in the long-wavelength limit, has been shown to exhibit fermion character. The fermion-like behavior of phonons discussed herein arises from the coupling between waves propagating along the positive and negative directions and from the projection of the wave function on the orthogonal space of direction of propagation. This fermion behavior is associated with constraints reminiscent of Pauli's Exclusion Principle. The present exclusion principle is a reflection of the exact nature of the eigenvectors of the wave equation projected on the space of directions. The eigenvectors of the phononic structures presented above are not simply plane waves but a superposition of plane waves traveling in opposite directions. These superposition states can be visualized as quasi-standing waves possessing a standing wave character and a traveling wave contribution. For a given wave vector, it is the fixed ratio of the amplitudes of the standing wave component and of the traveling wave components that constrains wave propagation and imposes an exclusion principle-type restriction. That is, an additional traveling wave cannot be superposed with the same wave vector without creating a new quasi-standing wave that is not solution of the wave equation.

One application of the properties of the wave function associated with its spinor character has been demonstrated. Encoding of information in the elastic waves that are supported by the phononic structures described above leads to the potential of operating simultaneously on the components of the spinor part of the wave functions. These operations act on the phase of states that can be written as superposition in a basis of pure forward or backward direction state and are called φ-bits. The present disclosure provides opportunities in operating on the states of multiple φ-bits. in a parallel fashion.

Finally, it is necessary to list some of the challenges that may be encountered in realizing experimentally a device that exhibits the fermion-like behavior described herein. These challenges are associated with (a) the difficulties in exciting rotational elastic waves, (b) difficulties in exciting only pure rotational modes without exciting other longitudinal or transverse modes, (c) loss resulting from the intrinsic viscoelastic properties of the media constituting the device, and (d) loss resulting from geometrical non-linearity associated with rotational degrees of freedom. 

The invention claimed is:
 1. A method for encoding information in a phononic structure, the method comprising: providing a phononic structure having a fermion-like behavior and supporting elastic waves, the phononic structure describable by a wave function, wherein the wave function has a spinor part, the spinor part having one or more components; and modifying the one or more components of the spinor part, wherein modifications of the components of the spinor part act on phase-bits, each phase-bit having a state comprising a superposition of waves propagating in a direction through the phononic structure; wherein the state of the phase-bits is associated with information encoded therein.
 2. The method of claim 1, wherein the waves only propagate in a pure forward direction and a pure backward direction.
 3. The method of claim 1, wherein the phononic structure supports rotational waves characterized by wave functions including a spinor part and a spatio-temporal part.
 4. The method of claim 1, wherein the phononic structure comprises springs, the springs comprising an elastic material.
 5. The method of claim 4, wherein the elastic material has a stiffness, the stiffness dependent on the magnitude of an externally applied electromagnetic wave.
 6. The method of claim 4, wherein the elastic material has a stiffness, the stiffness dependent on the magnitude of an external magnetic field.
 7. A system for encoding information in a phononic structure, the system comprising: a phononic structure having a fermion-like behavior and supporting elastic waves, the phononic structure describable by a wave function, wherein the wave function has a spinor part and the spinor part has one or more components; and one or more phase-bits having a state comprising a superposition of waves propagating in a direction through the phononic structure, wherein operating on the components of the spinor part of the wave function acts on the acts on the state of the phase-bits; wherein the state of the phase-bits is associated with information encoded therein.
 8. The system of claim 7, wherein the waves only propagate in a pure forward direction and a pure backward direction.
 9. The system of claim 7, wherein the phononic structure supports rotational waves characterized by wave functions including a spinor part and a spatio-temporal part.
 10. The system of claim 7, wherein the phononic structure comprises springs comprising an elastic material.
 11. The system of claim 10, wherein the elastic material exhibits photo-elastic effects.
 12. The system of claim 10, wherein the elastic material exhibits magneto-elastic effects.
 13. A system for encoding and operating on information in a phononic structure, the system comprising: a phononic structure having a fermion-like behavior and supporting elastic waves, the phononic structure describable by a wave function, wherein the wave function has a spinor part and the spinor part has one or more components; two or more phase-bits having a state comprising a superposition of waves propagating in a direction through the phononic structure, wherein operating on the components of the spinor part of the wave function acts on the state of the phase-bits, wherein the state of the phase-bits is associated with information encoded therein; and a device linking a magnitude of an external stimulus on a second phase-bit to a transmission coefficient of a first phase-bit, wherein two external stimuli are applied on the first phase-bit and the second phase-bit simultaneously by launching a wave into the phononic structure.
 14. The system of claim 13, wherein the waves only propagate in a pure forward direction and a pure backward direction.
 15. The system of claim 13, wherein the phononic structure supports rotational waves characterized by wave functions including a spinor part and a spatio-temporal part.
 16. The system of claim 13, wherein the phononic structure comprises springs comprising an elastic material.
 17. The system of claim 16, wherein the elastic material exhibits photo-elastic effects.
 18. The system of claim 16, wherein the elastic material exhibits magneto-elastic effects. 